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O ' Abstract 

J3 ' An SIRS epidemiological model for a vertically transmitted disease is discussed. 

We give a complete global analysis in terms of three explicit threshold parame- 
ters which respectively govern the existence and stability of an endemic proportion 
^^ . equilibrium, the increase of the total population and the growth of the infective 

CN ■ population. This paper generalizes the results of Busenberg and van den Driessche. 
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X ■ 1 Introduction 

In 1990, a complete global analysis of an SIRS model of disease transmission in a pop- 
ulation with varying size was given by Busenberg and van den Driessch 0]. In an SIRS 
epidemiological model, we divide the population into three groups, Susceptibles, Infec- 
tives and Removeds, and the problem is to examine the behavior of the size of each group 
when the time goes to infinity. They considered a disease with horizonal transmission, 
that is a disease which is transmitted by contact between an infective and a susceptible 
individual. We intend to generalize their results for a vertically transmitted disease, that 
is a disease which is also transmitted from infective parents to their newborns. The as- 
sumption of vertical transmission has two consequences. First it causes some newborns 



to die and forces us to assume that the birth rate differs from one group to another. (See 
the demographic assumption in |^, ^ [T^-) The second fact is that some newborns are 



infected, hence a group of newborns enter to the infective class [|l], ^ § . We also assume 
that a part of these infected newborns are known and removed after their birth. We shall 
show that the latter parameter can play an important role in the epidemic process. 

We consider an SIRS epidemiological model for a vertically transmitted disease. In 
our model, the incidence function is of proportionate mixing type introduced by Nold 



Tl| . Natural births and deaths are assumed to be proportional to the class numbers with 
different rates. We also assume that a proportion of the infected newborns are known and 
removed. We follow 0] to examine our model equation which is homogeneous of degree 
one. We consider the proportions system and show that this system has no periodic 
orbit in its feasibility region. This reduces our analysis to the discussion of existence and 
stability of rest points of a palnar system. The technique used here to show the uniquness 
of endemic equilibria is based on the Poincare index. This technique has no hard analysis 



and can be easily applied to other similar systems [|T3| , |1J]. The reader can verify that 
our results hold for similar SIRI systems as well [0]. 

We first in the next section state the model and a result concerning the non-existence 
of certain types of solutions . We consider the proportions system and prove that every 
solution in the feasibility region tends to a rest point of this system. In Section 3., we 
introduce a quadratic planar system with the same dynamics as the proportions system 
and then we discuss the existence and stability of rest points of this quadratic planar 
system. This gives a complete global analysis of the proportions system which is used to 
provide a global analysis of the original system is Section 4. 

2 The model 

In order to derive our model equations, we divide the population into three classes, the 
susceptible, the infective and the removed individuals with total numbers S, I and R 
respectively,. We set A^ = S + I + R which is the total size of the population. The 
following parameters are used in our model equations: 

bo: per capita birth rate of susceptible individuals. 



bi: per capita birth rate of infective individuals born uninfected, 

/3: per capita birth rate of infective individuals born infected, 

62: per capita birth rate of removed individuals, 

d: per capita disease free death rate. 

El. excess per capita death rate of infected individuals, 

€2: excess per capita death rate of removed individuals, 

a: per capita removal rate of infective individuals, 

7: per capita recovery rate of removed individuals, 

A: effective per capita contact rate of infective individuals. 

As mentioned before, we assume that the infected newborns enter the classes / and R 
of proportions (3i and /?2 respectively, hence (3 = Pi + (32- In this paper, all the above 
parameters, are positive, however some of them can also be zero. These hypotheses yield 
the following system of differential equation in M.'^ , where " ' " denotes the derivatives 
with respect to t, the time. 

^ S' = {bo- d)S + bj + {b2 + ^)R- \^-^ (2-1) 

I' = {f3i-d-ei-a)I + X— (2-2) 

^ R' = {(32 + a)I -{d + e2 + -f)R (2 - 3) 

IS 
where A—- is of the proportionate or random mixing type P, |lT|. By adding the above 

three equations, the total population equation is 

N = boS+{bi + f3- ei)I + (62 - e2)R - dN. 

SIR 
If we consider the proportions s = — , i = — , r = — , we get the following system of 
t^ f N N N 

equations. 

s' = bos + bii + (62 + l)r - bos'^ - (61 + /3 + A - 6i)is - (62 - £2)sr (2 - 1)' 

i' = (A -ei- a)i + (A - bo)is -{bi + (3- £1) ^ - (62 - £2) ir (2 - 2)' 

r' = {(32 + a)i - (£2 + l)r - b^sr - {bi + (3 - ei)ir - (f>2 - £2)r^ (2 - 3)' 
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This system is called proportions system and the feasibility region of this system is the 
triangle 

D = {{s,i,r)\s > 0,z >0,r > 0,s + z + r = 1}. 

If we set S = s + i + r then S' = (1 — S)(6oS + (61 +/? — £i)i + (62 — £2)'^)- Thus the plane 
S = 1 is invariant. Moreover, on the sides of D, we have: 



s' = bii + (62 + l)r > 0, 

i' = hence the line {i = 0} is invariant, 

r' = (/?2 + a)i. 



Therefore, D is positively invariant. On the invariant line {S = 1} fl {i = 0}, we have 

r' = -{e2 + 7)r - bor{l - r) - (62 - £2)^^ = -(^o + £2 + lY + {bo - &2 + £2)^^- 

It follows that this invariant line contains two rest points, the Disease-Free Equilibrium 
(1,0,0) and possibly another one which is outside of D. It is easy to see that the DFE 
attracts the side Dr\{i = 0}. Furthermore, our vector field is strictly inward on the other 

o 

sides of D. Thus D, the interior of D, is positively invariant too. The following theorem 
reduces our problem to the discussion of existence and stability of rest points in D. 



Theorem 2.1. The a;-limit set of any solution for the system (2 — 1)' — (2 — 3)' with 
initial point in D is a rest point in D. 



Proof. Since D is compact and positively invariant, the w-limit set of any solution with 

initial point in D is a compact nonempty invariant subset of D. Here we use the Poincare- 

Bendixon theorem and the terminology used to prove it [|12]. In the next section, we will 

see that this system has at most three rest points in D. Thus it satisfies the assumptions 

of Poincare-Bendixon theorem. We follow [^ |^ and define the vector field g = {gi,g2, gs) 

on D by 
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where /i, /2 and /a are the right hand side of (2 — 1)', (2 — 2)' and (2 — 3)' reduced to 

o 

functions of two variables by using ^ = 1 respectively. Clearly g.f = in D and after 
some computations 0], we get 

(curl g).{l, 1,1) = -^^ + — -— + — 

\s r SI sr"^ 

Since the DFE is the only invariant subset of dD (i.e. the boundary of D), the cj-limit 

o 

set must have some regular point in D if it is not a rest point. Let x be such a regular 
point and h be the first return map (Poincare map) defined on a tranversal at x. For a 
point y near x on the transversal, Let V be the region surrounded by the orbit F from y 
to hijj) and the segment between them. (This region is known as Bendixon Sack, Sec Fig 
2.1) Now by Stokes' theorem 

(curl^).(l,l,l)(ia= I g.fdt+ [ g{ty + {I -t)h{y)).{y - h{y))dt 
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Since g.f = and h{x) = x, the right hand side of the above equality tends to zero when 
y tends to x, but the left hand side tends to the integral over the region bounded by the 

o 

a;-limit set, this is a contradiction since (curl g)-{l, 1, 1) < in D- □ 




Fig. 2.1. The Bendixon Sack. 



3 The planar system 

using the relation s+i + r = l,we see that our system is essentially two dimensional. Thus 
we can eliminate one of the variable to arrive at the following quadratic planar system: 

i' = {X + Pi-bo-Ei- a)i + (bo + Si - X - h - (5)f + (feo + £2 - A - b2)ir (3 - 1) 
r' = {P2 + a)i - (bo + £2 + l)r + (&o + £1 - bi - (3)ir + {bo + £2- 62)^^ (3 - 2) 



The dynamics of the system (2 — 1)' — (2 — 3') in D is equivalent to the dynamics of this 
planar system in the positively invariant region Di = {{i,r) : z > 0, r > 0, i + r < 1}. 
This quadratic system has at most four rest points and since Di misses a rest point on 
i = 0, there are at most three rest points in Di. One of these rest points is the origin 
which comes from the DFE. The matrix of the linearization of the system (3 — 1), (3 — 2) 
at the origin is: 

\ + l3i — bo — El — a 

(32 + a -{bo + £2 + 1) 

with the eigenvalues X + Pi — bo — £1 — a and —(60 + £2 + t)- Now we define the first 

threshold parameter Ro = which governs the stability of the origin. 

bo + £1 + a 

Theorem 3.1. The origin is globally asymptotically stable in the feasibility region Di 
when -Ro ^ 1 and it is a saddle point when Ro > 1. 



Proof. By Theorem 2.1. it is enough to prove that ii Ro < 1, then the origin is the only 

o 

rest point in Di. If there exists a rest point in Di, we have i' = and i 7^ at this point. 
Thus it belongs to the line 

{X + Pi-bo-ei-a) + {bo + ei- X-bi- p)i + (bo + £2 - X - 62)^ = 0. (3 - 3) 

Since s' = at this point, from (2 — 1)' we obtain 

boS + bii + (62 + 7)^ - bos'^ -{bi + (3 + X- ei)is - (^2 - £2)sr = 0. 

and by using the relation s + i + r = 1, we can write 

bii + (62 + l)r + {bo + £i-bi- (3- X)is + {bo + £2 - &2)sr = 0. 

Multiplying (3-3) by (— s) and adding it to the above expression, we get the following 
equality 

bii + (62 + l)r + Xsr + {bo + £i + a- X- (3i)s = 0. 

But the left hand side is positive when Ro ^ I and this is a contradiction. D 

When i?o > 1, the origin is a saddle point and it does not attract any point of 
Di — {i = 0}. Thus the orbits with initial point in Di — {i = 0} must be attracted by 



some rest points in Di by Theorem 2.1. These rest points belong to the hne (3-3) and 
the conic section 

r' = {P2 + a)i - (bo + 62 + 7)^ + (&o + £1 - ^i - /3)ir + (60 + £2 - ^2)^^ = 0. 

o 

It follows that there are at most two rest points in Di- Notice that a nondegenerate rest 
point of the planar system is obtained by a transverse intersection of the line (3-3) and 
the above conic section. The following lemma has two immediate consequences which will 
be very helpful. 

Lemma 3.2. The trace of the linearization of the system (3-1), (3-2) at a rest point in 

o 

Z^i is negative. 

Proof: We compute the trace at a rest point. 
di' 



{\ + Pi-bo-Ei-a) + 2(60 + 81- X-bi- f3)i + {bo + e2- X- 62)^ 



^ = -(60 + £2 + 7) + (^0 + ei-bi- p)i + 2(60 + €2- b2)r. 
or 

From i' = and r' = at a rest point, we can write 

0% uv t 

-— = {bo + ei- X-bi- (3)i and t^ = (^0 + £2 - ^2)^ - (/32 + «)-• 
m or r 

Since s' = bii + (62 + '^)r + {pQ + £1 — X — bi — [3)si + (60 + £2 — b2)sr, it follows that 

di' Or' 
(bo + £1 — X — bi — [3)si + (60 + £1 — A — 61 — I3)sr < and hence — - + — — < 0. D 

oi or 



Corollary 3.3. The system (3-1), (3-2) has no source point in Di 

o 

Corollary 3.4. Every nondegenerate rest point in Di is hyperbolic. 

o 

Theorem 3.5. If Rq > 1, then there exists a unique rest point (i*,r*) in Di which is 
hyperbolic and attracts Di — {i = 0} . 



Proof. When i?o > 1, the origin is a saddle point with the unstable eigenvector 

X + Pi - bo - ei - a + bo + 62 + 'y 
(32 + a 

Since -Rq > 1, we have X + (ii — bo — Si — a > Q and hence this vector belongs to the first 

o 

quadrant of the plane (i,r). Since Di is positively invariant it follows that a branch of 

o 

the unstable manifold of the origin lies in Di- (See Figure 3.1.) 




Figure 3.1. Local behavior of planar system near the origin when i?o > 1 

This helps us to find a piece-wise smooth Jordan curve C on which our vector field X is 
either tangent or inward. (See Figure 3.2.) The Poincare index of such a Jordan curve is 



1. (See [jT3l, Lemma 5.1.) 




Figure 3.2. The Jordan curve C 

o 

We choose this Jordan curve so that it contains all rest points in Di- If there are two 

o 

rest points in Di^ they are obtain by a transverse intersection of the line (3-3) and the 
conic section r' = and hence both are nondegenerate. Thus their Poincare index must 

o 

be ±1 which contradicts Ix{,C) = 1. Therefore there is a unique rest point in Di which 
attracts Di~{i = 0} by Theorem 2.1. It remains to prove that this rest point is hyperbolic. 
Suppose the contratry, then it must be nondegenerate by Corollary 3.4. Thus it is obtained 
by a tangent (non-transverse) intersection of the line (3-3) and the conic section r' = 0. 
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With a slight perturbation in 7, we will have either two transverse intersection in Di or 
nothing. (7 appears only in the coefficient of r in r' = and does not appear in (3-3) and 

o 

Rq). This is a contradiction with the uniqueness of the rest point in Di proven above. D 
Theorem 3.1. and Theorem 3.5. provide a complete global analysis of the planar 

o o 

system (3-1), (3-2) in Di- Since the dynamics of this system in Di is equivalent to the 
dynamics of proportions system (2 — 1)' — (2 — 3)' in D, we have proved the following 
result which gives a complete global analysis of the proportions system in the feasibility 
region D. 

Theorem 3.6. Consider the proportions system (2 — 1)' — (2 — 3)'. 

(i) If Rq < 1, then the disease free equilibrium proportions (1,0,0) is globally asymp- 
totically stable in D. 

(ii) If -Rq > 1) then there is a unique rest point (s*, i* , r*) which is globally asymptotically 
stable in D — {i = 0}. 



4 Analysis of the model equations 

Consider the original model equation (2-l)-(2-3) and recall that the population equation is 
N' = boS +{bi+/3- ei)I + (62 - e2)R - dN. Thus 

N' 

— = bos + {bi + f3- ei)i + (^2 - e2)r - d. (4 - 1) 

N' 
If Rq < 1, then (s,i,r) — > (1,0,0) by Theorem 3.6., hence > bg — d. Moreover if 

o 

i?o > 1 and / > 0, then (s,i,r) — > (s*,i*,r*), i.e. the unique rest point in Di, and 

N' 

— ^ 60S* + (61 + /? - e)i* + (62 - e2)r* - d. 

We define the second threshold parameter which governs the total population as follows. 

I if/2o<l, 

^1 - < bos* + (61 + (3)t* + b2r* .^ ^ ^ 
— — II Ro > 1- 

a + Ell* + e2r* 

9 



Since J = is invariant with a linear equation, we may assume that / > 0. Now from 

(2-2), we write 

_ = (/3i -rf-£i -a) + As. (4-2) 

If i?o < 1, then — — > Pi — d — El — a + X and if Rq > 1, — — > (3i — d — ci — a + As*. 
So we define the third threshold parameter which governs the total number of infective 
individuals. 

if ito < 1, 



T> — ) d + ei + a 

d + ei + a 

Notice that d does not appear in the proportions system and hence (s*, i*,r*) is indepen- 
dent of d. The following results provide a rather complete global analysis of the model 
equations (2-l)-(2-3). 



Lemma 4.1. If I(t) < M for every t > to, then R{t) — > and if I{t) — > cxd, then 
R{t) — > oo. 

Proof. From (2-3), we have R'{t) = {P2 + a)I{t) - {d + 62 + -f)R{t). If /(t) < M 
for t > to, then R'{t) < M{P2 + a) — {d + 62 + "y)R(t) and by Granvell's inequahty 
H, R{t) < M{P2 + a) e-('^+"^+T') (*-*") which follows that R{t) — > 0. Now suppose that 



I{t) — > 00. Then by (2-3), we have R'{t) + {d + e2 + l)R{t) = (/?2 + a)/(t) which implies 

d 
that — (i?(t)e('^+=2+^)*) = (/52 + a)/(t)e('^+"2+7)t. por every M e M+, there is a to G M such 

dt 
that I{t) > M for t > to. Thus 

^(^)g(rf+e2+7)i _ i?(to)e('='+"2+7)to = f {[32 + a)/(t)e(''+"2+^)* > M{p2 + a) f e^'^+^^+^^^dt. 

Jto -/to 

^ uit) > (ij(i„) - M(A+j^) ,-(.«..,)(.-.,.. + £Wi±i!) . 

^^ V d + £2 + 7/ ^ + £2 + 7 

Since e~*^'^+^2+7)(*-*o) croes to zero as t — > 00, R(t) > —, for large values of t. 

^ ' ^^ 2((t + £2+7) 

It means that i?(t) — > cxo. D 



Theorem 4.2. (i) If i?i > 1, then N{t) — > 00 and if i?i < 1, then N{t) — ^ 0. 
(ii) If i?2 > 1, then {I{t), R{t)) — > (00, 00) and if R2 < 1, then {I{t),R{t)) — > (0, 0). 

10 



Proof: First suppose that -Ro ^ 1 which imphes that {s,i,r) — > (1,0,0) by Theorem 
3.6. Thus — — > b- dhy (4-1) and - — > l3i - d - ei - a + X hy (4-2). If Ri < 1, 
then b- d <0, hence N{t) — > 0. Similarly if Ri > 1, then b - d > 0, hence N{t) — > 0. 
Furthermore If i?2 < 1, then Bi — d ~ Si — a + X < 0, hence I{t) — > and by the above 
lemma R{t) — > 0. Similarly if i?2 > 1, then j3i — d — Si — a + X > 0, hence I{t) — > oo 
and by the above lemma R{t) — > oo. 

Now suppose that Rq > 1. We assumed that I(t) > 0, hence z > and (s,z,r) — > 

N' 
(s*, i*,r*) by Theorem 3.6. Thus — — > bos* + (6i + /? - ei)i* + (62 - £2)r* - d which is 

positive if Ri > 1, hence N{t) — > 00 and negative if Ri < 1, hence N{t) — > 0. Moreover 

r 

— — ^ {(3i — d — El — a) + Xs* which is positive if /?2 > 1, hence I{t) — ^ 00 and then 
R{t) — > 00 by the above lemma. Similarly (3i — dsi — a + Xs* < if i?2 < 1 and then 
(/(t),i?(t))— ^(0,0). D 

We summarize our results in the following table which is the same as Table 1. in 0|. 



Ro 


Ri 


R2 


N — > 


{s,i,r) — > 


{S,I,R)-^ 


< 1 


< 1 


< 1° 





(1,0,0) 


(0,0,0) 


> 1 


< 1 


< 1° 





{s*,i*,r*) 


(0,0,0) 


< I 


> 1 


< 1 


00 


(1,0,0) 


(00,0,0) 


< I 


> 1 


> 1 


00 


(1,0,0) 


(00, 00, 00) 


> 1 


> 1 


> 1" 


00 


(s*,z*,r) 


(00, 00, 00) 



" Given Rq and Ri, this condition is automatically satisfied. 

Conclusion Remarks. Here we will make some comparison between the result in 0] 
and the above results. 

1. If we set bo = bi = b2 = b and /5i = /32 = /? = in the system (2-l)-(2-3), we obtain 
the system (2-1)- (2-3) in 0]. Moreover by these assumptions, we get Rq = , 



R, 



b 


ifi?o<l, 
if Ro>l 


d + Sii* + e2r* 
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,^ ^ ^ if Ro < 1, 

— ^^^ ifi?o>l. 

a + £i + a 

which are the same threshold parameters as in [Q. 

2. Comparing their threshold parameters with our ones, we see that the effect of 
60, hi and 62 appears more clearly instead of h. For example when i?o > 1, our Ri is 

°^ t ^t I "*", ^'' , but they obtain Ri = -r, — ^, — - in which the effects of hn, 61 and ho are 
hidden in 6. 

3. There are two vertical transmission parameter, (3i and /^i, in our model. The effect 
of (3i in Rq and R2 is crucial. Since /3i + /?2 = /5, we can decrease /3i by increasing j32. 
Therefore removing more infected newborns causes more safe situation. 

Acknowledgments. The author would like to thank Institute for Studies in Theoretical 
Physics and Mathematics for supporting this research. 
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